LETTER TO THE EDITOR 



43 



X 



A new description of spin tunneling in magnetic 
molecules 



D. Galetti 



O 

Institute de Ffsica Teorica, Universidade Estadual Paulista, Rua Pamplona 145, 

01405-900, Sao Paulo, SP, Brazil 
5-H ■ 
Oh' E-mail: ealetti@ift .unesp.br 

< 

l/~J ■ Abstract. A new approach is used that allows to describe the magnetic molecules 

main properties in a direct and simple way. Results obtained for the Fes cluster show 

T 1 

^. ' good agreement with the experimental data. 

^j- ' PACS numbers: 03.65. Ca, 75.45+j, 75.60Jp 

o" 

O 



Submitted to: J. Phys. C: Solid State Phys. 



In recent years, with the experimental advances in the measurement of magnetic 
molecular clusters properties, it has emerged a new frontier in this area. One of the 
interesting aspects that came out of these studies is the possibility of measuring spin 



tunneling in mesoscopic systems what corresponds, in a standard quantum description, 
to the tunneling of the collective degree of freedom corresponding to the magnetization 
direction through a potential barrier separating two minima of an effective potential 
associated with the spatial orientation [TJ El El HI El El 13 IB] • 

From the theoretical point of view, spin tunneling has been treated mainly by the 
use of a WKB method adapted to spin systems [H| HO], by using Feyman's path integral 
treatment of quantum mechanics |1H IT2]. and also by using su(2) coherent states J3] in 
order to establish a correspondence between the spectrum of the spin system with the 
energy levels of a particle moving in an effective potential[14|. 

In the present letter we intend to show that still another approach may be used 
for describing spin tunneling - in angle representation - in such a way that analytic 
expressions for the characteristic parameters of the magnetic molecular clusters can 
be obtained; in particular, the results for the spectrum and energy barrier heights are 
directly obtained and the spin tunneling process can be easily interpreted. 

The starting point of the present approach is the introduction of a quantum 
phenomenological Hamiltonian describing the spin system, written in terms of angular 
momentum operators obeying the standard commutation relations, and that reflects 
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the internal symmetries of the system. It may also contain terms taking into account 
external applied magnetic fields. The degree of freedom that undergoes tunneling is 
considered a particular collective manifestation of the system, and it is assumed to be 
the only relevant one. At the same time, the temperature of the system is assumed 
so conveniently low that possible related termally assisted processes are not taken into 
account so that only quantum effects are considered. For instance, the general quantum 
Hamiltonian - sometimes also called giant spin model - 

H = AJ z -DJ 2 z +^{Jl + J 2 _) (1) 

can be used to study some systems of interest. In particular, this Hamiltonian can 
describe the octanuclear iron cluster JT3], Feg, in the presence of an external magnetic 
field along the z axis if we take A = gfiBH\\] by taking A = no external field is 
considered. This spin system has a j = S = 10 ground state and a suggested pure 
quantum spin tunneling below 0.35 K\ furthermore D/ks = 0.275 K, E/kg = 0.046 K 
[UJIH], and ks is Boltzmann's constant respectively. On the other hand, from a pure 
algebraic model point of view, it is interesting to see that the Lipkin quasi-spin model [To] 
of wide use in many-body physics is also obtained by just considering D = 0. At the 
same time that the Hamiltonians in both cases are similar, and in this form we can 
compare their results, our interest in this model also resides in the fact that it stands 
for a valuable testing ground for checking the validity of approximations in treating 
collective degrees of freedom. 

In what follows we will show how we can discuss the spin tunneling process in the 
Fes cluster by the use of a new Hamiltonian that is an approximate version of Eq. ((TJ . 
This new Hamiltonian is obtained through a series of transformations performed on 
the matrix generated by calculating the expectation values of Eq. (£Q) with the su(2) 
coherent states \j,z), where z is a complex variable and j characterizes the angular 
momentum state rmiltiplet[Kl|. The Hamiltonian is the overcomplete spin coherent 
states representation then given by 

(j,z'\H\j,z) = K{z',z), (2) 

the also known generator coordinate energy kernel [T7j. embodies the quantum 
information related to the system we want to study. The procedure of extracting a 
new Hamiltonian - written now in terms of an angle variable - from Eq. (J2J) has been 
already shown elsewhere and will not be repeated here (THJ [THl 1201- In fact, as it was 
proved there, the variational generator coordinate method can in this case be used to 
rewrite the Hamiltonian from which we start in an exact and discrete representation 
which can then be conveniently treated in order to give an approximate Hamiltonian in 
the angle representation. It is important to point that the Hamiltonian we obtain is in 
fact an approximate one, but we have also shown, by studying the Lipkin model, that 
it is already a reliable Hamiltonian for spin systems with j = S > 5, as it is the case for 
the Fes cluster, when S = 10, as mentioned before. Furthermore, since this approach 
is based on quantum grounds from the beginning, it does not need to go through any 
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quantization process. Also, it is not necessary to convert the discrete spin system into 
a continuous one as it is usually done|21j. at the same time that the quantum character 
of the angle- angular momentum pair is properly taken into account. 

The Hamiltonian in the angle representation associated with the Fe$ cluster in the 
absence of an external magnetic field is then written explicitely as 

fl '« = i« + KW - (3) 

where 

V(<j)) = -(D-E)S{S + 1) cos 2 <f)-ES{S + l) (4) 

is the potential energy, whereas the effective "mass" is given by 

M ^ = 2(D-E)L^ + 4E ' (5) 

with — 7r < <p < 7T. It is important to observe that the effective "mass" is not constant 
over the angle domain and that it plays an essential role in this description. Figure 
1 shows the potential function as well as the effective "mass". The minima of both 
functions occur at (f> = 0, tt while the maxima occur at <f> = |, %-(= — f ), as expected. 

The eigenvalues and eigenfunctions associated with the Hamiltonian Q can then 
be directly obtained by numerically solving the Schrodinger equation 

H(<f ) )ij k ( ( j ) ) = S k ^ k (<j ) ) (6) 

by means of a Fourier analysis. The ground state energy thus obtained, £ gs ~ 
—27.6447 K, agrees quite well with the result obtained by numerically diagonalizing 
the exact phenomenological Hamiltonian Eq. (£Q) within the \j = S, m) state basis 
(hereafter these results will be called the reference values), specifically, the deviation 
from the reference value is of the order of 0.5%. Figure 2 depicts the wave functions 
associated with the lowest pair of energy eigenstates. In what concerns the energy 
splitting of those states, we obtain a result that is of the same order of magnitude as 
the reference value; in fact, the deviation is of the order of 17.5%, being that our result 
is smaller than the reference value. 

Now, we see that the top of the potential barrier is explicitly given by 

V max (<j>) = -ES(S + l) (7) 

so that 

h b = -ES(S + l)-E gs (8) 

measures the ground state energy barrier. Using our result for the ground state we get 
hf, ~ 22.58 K, which is only 1.7% higher than the experimental result, namely 22.2 K, 
presented in [6J. An approximate analytic expression for the ground state energy barrier 
can be obtained from another perspective. To this end, we first take into account that 

£mm = V ((f) = 0) = -DS (S + 1) = -30.25 K, (9) 
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and we also assume that the ground state energy is approximately given by S gs 
£ m in + w/2. We then perform a harmonic approximation at the potential minimum, 

d 2 V (0) . 



M(0)a; 2 U (Ess) 

UlfJ 

from which we obtain 

2yJ{D-E)[ES{S 



\£gs\)- 



Taking advantage of this dependence on \S gs \ we obtain the analytic expressions 

D F 



\S gs \ = DS(S+l] 



and 

h b = {D-E)S(S 



E 



D + E 
D-E 



S(S + 1) 



27.52 K, 



1) 



D-E 



1 + 4- 



D + E 



S{S + l) 



22.46 K 



(10) 



:n^ 



(12) 



(13) 



D-E 

for the ground state energy and barrier respectively. The deviation in the ground state 
energy is then 0.08% while for the barrier height it is 1.17%. Even if we consider the 
crude approximation (D + E) / (D — E) ~ 1.0 we get 

\E gs \ ~DS 2 + ES = 27.9QK, (14) 

and 

h b ~(D- E) S 2 = 22.90 K, (15) 

respectively; the value for the barrier height still is in good agreement with the 
experimental result. 

Now, if an external magnetic field paralel to the z axis is applied the potential 
function reads 

V((f)) = -(D-E)S(S+ 1) cos 2 <f>-y/S(S+l)gfj, B H\\ cos 0-ES (S + 1) , (16) 

while the expression for the effective "mass" is given by 

1 
M ^ = 2 (D-E) cos 2 + ^#11 cos0 + IE' ^ 

It is immediate to see that the presence of the magnetic field does not change the position 
of the potential minima, = 0, ir, but it introduces a shift in energy at these points so 
that the difference in the height between the two minima will be (H\\ > 0) 

V 



n 



V (0 = 0) = 2y/S(S+l) gf x B H t 



11 



This means that for some particular values of H\\ there will occur a degeneracy in the 
energy spectrum such that the state with, for instance, m = n will match its energy with 
that of the state with m! = —n + k, as it is indeed expected. The particular value of the 
magnetic field increment Ho that leads to the matching of the energy levels, therefore 
given rise to the appearance of the degeneracies, H\\ = k H , can be obtained from an 
analysis of the effective "mass" expression. Realizing that the presence of zeroes in the 
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function 1(0) = 1/M (0) (infinities of the effective "mass") indicates that tunneling 
cannot occur, we look for the expression for the strength of H\\ beyond which tunneling 
will not take place. A direct calculation shows that the limit is given by 

fl-lim = ^SHb ^ 2E p _ E j _ 4 32 T ^ 

so that 

H}r 2hu , 

H = — L = —W2E (D-E)~ 0.216 T (20) 

26 gfx B 

for g = 2. This result is in good agreement with the experimental value H Q = 

0.22 T[6, 8J. It can be immediately seen that for this value of the external paralel 

magnetic field the original minimum at = it and the maxima at = it/2 and = 37r/2 

have turned into a single maximum of the potential function, while the only surviving 

minimum is the one at = (= 2ir). This means that, in this particular situation, 

there is only one direction the spin can be directed at. 

In conclusion we have proposed an angle-based description of the spin tunneling 
in magnetic molecules that can account for the basic results governing this kind of 
phenomenon. The results obtained in the case of the Fes cluster agree quite well with 
the experimental data and the interpretations follow in a direct and simple way. 
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Figure Captions 

Figure 1: The potential and effective "mass" functions associated with the Fes 
cluster with parameters D/ks = 0.275 K and E/ks = 0.046 K. The interval is taken 
as —tt/2 < <p < 3tt/2 for clarity. 

Figure 2: Fe§ cluster ground state wave function if) ((f)) and first excited state wave 
function ipi (0). 
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